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Jet masses for high p T QCD jets at all orders. 
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The study of the shape and sub-structure of high pj jets produced in hadron collisions is becom- 
ing an increasingly important component of LHC phenomenology in the context of new particle 
discoveries. We study here the state of the art for masses of QCD jets in perturbation theory to 
all orders, up to next-to-leading logarithmic accuracy, with inclusion of non-global logarithms. 
We also study the precise role of the jet algorithm used to construct the jets and its impact on re- 
summed predictions. Such analytical predictions, where available, offer an alternative to leading 
logarithmic predictions from Monte Carlo event generators and could, after fixed-order matching 
and accounting for non-perturbative effects, be used for direct comparison to early LHC data. 
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1. Introduction 

Studies of jet-shape variables and in particular jet masses have long been a standard part of 
the program of testing perturbative QCD. The infrared and collinear safety of such observables has 
rendered them calculable by perturbative methods both in fixed-order perturbation theory as well 
as amenable to all-order resummed treatments in special regions of phase space enhanced by large 
logarithms [jlj]. 

In the context of the LHC, jet-shape and sub-structure techniques [||] are also rapidly gaining 
impetus for example as a means of discovering new particles that are highly boosted and hence 
whose decay products may end up in a single jet. As examples of recent sub-structure methods 
in this context one can refer to the development of jet filtering |^] and pruning techniques [Q] 
amongst a variety of tools being developed to reduce QCD jet backgrounds as well as to enable 
better reconstruction of mass-peaks for new heavy particles. 

In this paper we concentrate on what is known about jet shapes and in particular jet masses 
for QCD jets in all orders resummed perturbation theory. At the LHC one may be interested 
for instance in the mass or shape of one or more jets produced in high pj jet events of a given 
jet multiplicity. Recently it was suggested to study such observables with the aim of resumming 
logarithms in the shape variables for certain jets in multi-jet events while leaving the shape of 
other jets unmeasured []5| |0. Resummed predictions were provided aiming at next-to-leading or 
single-logarithmic accuracy in the shape variable distributions, within the context of soft-collinear 
effective theory (SCET). However as was emphasised in Ref. [|7|] single logarithmic accuracy for 
such observables necessitates a study and inclusion of non-global logarithms [|8j, ^|. Below we shall 
examine this issue and comment on the size of the non-global effect as well as looking at the role 
of the jet algorithm which is also significant to single-logarithmic accuracy. 

2. Individual jet masses in multijet events and non-global logarithms 

We can consider either the jet-mass in single-jet inclusive studies or focus on events of fixed 
multiplicity where we may pick a subset of all jets for study, as was the case in Refs. [§, §]. 

To illustrate our points here we examine a simple situation involving high pj dijet production 
where we measure the mass of one of the jets while leaving the other unmeasured. We can take 
the jets to be produced at zero rapidity wrt the beam without changing our conclusions. We shall 
also consider the limit of narrow well-separated jets in the sense that the jet radius R can be taken 
to be small compared to the interjet separation Ajj = 1 — cosfyy. We shall thus neglect terms of 
relative order /? 2 /A, ; . In this limit we find that owing to QCD coherence a rather simple picture of 
jet evolution emerges njn. In particular initial state radiation only contributes terms to the jet mass 
that vanish in the R — > limit and hence can be ignored in our approximation. We shall carry out 
fixed-order calculations below, in the soft approximation, both at LO and NLO accuracy and from 
these calculations infer the form of our resummed results. 

3. LO estimate in soft approximation 

A simple eikonal leading-order calculation for the cross-section for normalised jet mass 4M 2 /Q 2 
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to be below some value p [^] produces 

C F a s fQl 2 dco, f„ooR 2 \ C f a Sl 2 ^ 2 ~/„2 \ 
7T ip e /2^ (o \ Qp J 2% p v F; ' 

where Q is the hard scale of the process, in this case the jet p T . The above result holds no surprises 
- it is the familiar double-log estimate for jet-mass, with R being the jet radius. Such double loga- 
rithmic terms exponentiate and produce the Sudakov peak in the mass distribution. The calculation 
can be easily extended to account for hard emission collinear to the jet, which is a relevant source 
of single logarithms. However for the jet mass distribution at single-log accuracy there is more to 
the story than exponentiation of a single gluon. We report on this in the next section but note that 
at the level of terms in the conventional Sudakov exponent there is no involvement of the details of 
the algorithm except the radius R - all algorithms are identical at this level. 



4. NLO estimates and non- global logarithms 

Beyond leading order there is the issue of the soft gluon emission pattern as well as its interplay 
with the jet algorithm both of which lead to the appearance of non-trivial effects. A simplified 
emission pattern where soft gluons are coupled directly to the hard emitting ensemble is sufficient 
at the Sudakov level. However as was shown in Refs. |g[ ||] for non-global observables that are 
affected by soft radiation in a delimited phase space region such as the interior of a jet the pattern 
of real-virtual cancellations that ensure the exponentiation of a single-gluon is spoiled. One ends up 
needing to consider multiple soft gluon emission from an arbitrarily complex ensemble involving 
not just the hard partons but also all subsequent soft emissions until one hits the veto scale (the jet 
mass). This complication starts at the two-gluon level with the soft correlated two gluon emission 
term which for the quark dijets we assume here, has a CfCa colour factor. 

Moreover from the two gluon level onwards the details of the jet algorithm also become im- 
portant due to soft gluon clustering effects. One algorithm where these effects can be ignored is the 
anti-k, algorithm [10] in which soft gluons cluster to the hard jets long before their self-clustering 
can take place, leading to circular jets. In this algorithm it proves possible, at least in the large N c 
limit, to address relatively easily the complex dynamics leading to non-global logarithms. Carrying 
out the calculation for the emission of two soft gluons and focussing on the CfCa piece we obtain 

The above equation needs some explanation. Firstly we have introduced a parameter Eq analo- 
gous to the parameter A in Refs. |6|], which corresponds to the maximum energy flowing outside 
the hard high pj jets. Limiting Eq thus corresponds to restricting the hard jet multiplicity while for 
the inclusive jet mass case Eq should be understood to be of the order of the hard scale i.e typical 
jet pr- Here we shall focus on the inclusive jet mass case Eq ~ pj while in Ref. [Q] we also dealt 
with resummation of logarithms in pj /Eq. We note that the above equation is correct with neglect 
of additional terms suppressed by powers of R which are beyond our approximation. With this ap- 
proximation in place the coefficient 7T 2 /3 emerges and is the same coefficient as the corresponding 
one that appears for the hemisphere jet masses. This fact is not a coincidence but a result of the 
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Figure 1: Jet mass distribution with (blue curve) and without non-global logarithms (red curve) with Eq = 60 
GeV and R = 0.4. 

dominance of the collinear singularity that appears due to l/iki.bi) terms in the squared matrix 
element for the emission of soft gluons k\ and lc2- The non-global logarithms come essentially 
from the edge of the jet and have an additional analytic dependence on the jet radius R which can 
be ignored in a small R approximation. 

In the case of the hemisphere jet-mass studied in Ref. []|] non-global logs come from a large- 
angle boundary i.e that between the hemispheres. In the present case this boundary is replaced by 
the circular jet boundary but the coefficients and resummation of non-global terms will be the same 
as for the hemisphere case up to terms that vanish with R. The resummed result for the non-global 
piece, obtained numerically and in the large N c limit, can then be essentially taken from Ref. []8|] 
and reads 

K 2 2 ( \ + (at) 2 ' 



with t being the single-logarithmic evolution variable t = ^ f e - L ^ot s { x Q) and where L = In (2E R 2 ) / (Qp), 
a = 0.85Ca, b = 0.86Ca, c = 1.33. The non-global factor S multiplies the result obtained by ex- 



ponentiating the one-gluon result ( |3.1[ ) after taking account of hard collinear emissions and the 
running coupling. 

The role of non-global logarithms in the jet mass distribution can be seen from Fig. [|. The 
figure is for Q = 500 GeV which translates here to a 250 GeV jet pj with R = 0.4 and Eq = 60 GeV. 
There is roughly a twenty percent effect on peak height at these values of the parameters and in fact 
we find no significant variation in this number when Eq is lowered due to the increasing importance 
of non-global logs in pj /Eq as reported in Ref. The general situation where one measures the 
shapes of one or more jets in multi-jet production at high transverse momenta can be addressed by 
modifying the resummed result based on naive single-gluon exponentiation with a product of form 
factors from each jet, each of which has the form of S(t) reported above. The above form is correct 
up to terms that vanish as powers of R 2 /A/ ; - with R the jet radius and Ay the interjet separation. 
There is a clear physical reason for the emergence of the simple results reported here. In the small 
R limit each jet evolves independently of others due to the dominance of collinear emissions - the 
interesting physics (NLL contribution to jet masses) comes separately from the boundary of each 
jet which falls, at small R, within the collinear regime. 

The results derived in this article are valid only in the anti-&, algorithm where soft gluons 
cluster to the hard jets independently of one another resulting in circular jets. This situation is 
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modified when one moves to describing jets in other algorithms, which we shall discuss below. 



5. Other jet algorithms 

Let us now comment on the situation in jet algorithms other than the anti-& ; . Focussing on 
sequential recombination algorithms such as the Cambridge- Aachen (C-A) [ Jfl] , [12| ] and the k t al- 
gorithm Jl3| , p4| ] we encounter beyond leading double logarithms a much more complex situation. 



Appleby and Seymour [15] were the first to point out and numerically study the role of k t clustering 
on non-global logarithms (for the case of gaps between jets). Their conclusions were that the effect 
of k t clustering reduced the non-global component because a harder gluon k\ is capable of cluster- 
ing a softer emission &2 in configurations that may otherwise have generated significant non-global 
effects. Similar conclusions will apply for the C-A algorithm. It was subsequently found that clus- 
tering also had an impact at NLL accuracy on real- virtual cancellations in the independent emission 
(as opposed to correlated emission) or global term [16]. While in Refs. [ 16|, 17] the clustering gave 



terms that vanished with R, in the present case owing to the collinear singularities relevant for jets 
at small R, the effect is independent of R. We find that for a quark jet the contribution of clustering 
to the p distribution starts at ff(a s ) 2 and the first such term reads 

A E cluster = _Q ?28C 2 2 I ln I (5 1} 

dp 2 b \2%) p p 

which is a relevant term for resummations aiming at single-log accuracy in the mass cross-section 
£(p) and arises purely due to the clustering inherent in the algorithm. Likewise similar effects will 
be present for the C-A case. Single logarithmic terms generated by clustering were shown to be 



resummable for the gaps between jets case [ |17p but a similar calculation has not yet been carried 
out for the present case of jet masses. 



6. Conclusions 

In conclusion we point out that NLL resummed predictions can be made for quantities in- 
volving the shapes of jets produced in complex multi-jet events and using the anti-&, jet definition, 
which include non-global logarithms in the large N c limit. Up to correction terms varying as pow- 
ers of the jet radius there is an independent non-global factor arising from the boundary of each 
measured jet. The product of such factors multiplies the exponentiated "single-gluon" result which 
is relatively straightforward to obtain. We do not expect neglected terms varying as powers of R or 
those suppressed as 1 /N 2 to make a visible difference to our predictions even for R as large as 0.7 
Phenomenological investigations can thus be carried out using the analytical results after matching 
to fixed-order to obtain the high-mass tail and assessing the role of non-perturbative effects. 



References 

[1] S. Catani, L. Trentadue, G. Turnock and B. R. Webber, "Resummation of large logarithms in e + e~ 
event shape distributions," Nucl. Phys. B 407 (1993) 3. 

[2] M. H. Seymour, "Searches for new particles using cone and cluster jet algorithms: A Comparative 
study," Z. Phys. C 62 (1994) 127. 



5 



QCD jet masses 



Mrinal Dasgupta 



[3] J. M. Butterworth, A. R. Davison, M. Rubin and G. P. Salam, "Jet substructure as a new Higgs search 
channel at the LHC," Phys. Rev. Lett. 100 (2008) 242001 [arXiv:0802.2470 [hep-ph]]. 

[4] S. D. Ellis, C. K. Vermilion and J. R. Walsh, "Techniques for improved heavy particle searches with 
jet substructure," Phys. Rev. D 80 (2009) 051501 [arXiv:0903.5081 [hep-ph]]. 

[5] S. D. Ellis, A. Hornig, C. Lee, C. K. Vermilion and J. R. Walsh, "Consistent Factorization of Jet 
Observables in Exclusive Multijet Cross-Sections," arXiv:0912.0262 [hep-ph]. 

[6] S. D. Ellis, A. Hornig, C. Lee, C. K. Vermilion and J. R. Walsh, "Jet Shapes and Jet Algorithms in 
SCET," arXiv: 1001. 0014 [hep-ph]. 

[7] A. Banfi, M. Dasgupta, K. Khelifa-Kerfa and S. Marzani, "High-pr jet shapes, non-global 
logarithms and jet algorithms," arXiv: 1004.3483 [hep-ph]. 

[8] M. Dasgupta and G. P. Salam, "Resummation of non-global QCD observables," Phys. Lett. B 512 
(2001) 323 [arXiv:hep-ph/0104277]. 

[9] M. Dasgupta and G. P. Salam, "Accounting for coherence in interjet E(t) flow: A case study," JHEP 
0203 (2002) 017 [arXiv:hep-ph/0203009]. 

[10] M. Cacciari, G. P. Salam and G. Soyez, "The anti-fc, jet clustering algorithm," JHEP 0804 (2008) 063 
[arXiv:0802.1 189 [hep-ph]]. 

[11] Y. L. Dokshitzer, G.D. Leder, S. Moretti and B.R. Webber, "Better jet clustering algorithms", JHEP 
9708 (1997) 01 [arXiv:hep-ph/9707323] 

[12] M. Wobisch and T. Wengler, "Hadronization corrections to jet cross sections in deep-inelastic 
scattering," arXiv:hep-ph/9907280. 

[13] S. Catani, Y. L. Dokshitzer, M. H. Seymour and B. R. Webber, "Longitudinally invariant K t clustering 
algorithms for hadron hadron collisions," Nucl. Phys. B 406 (1993) 187. 

[14] S. D. Ellis and D. E. Soper, "Successive combination jet algorithm for hadron collisions," Phys. Rev. 
D 48 (1993) 3160 [arXiv:hep-ph/9305266]. 

[15] R. B. Appleby and M. H. Seymour, "Non-global logarithms in inter-jet energy flow with k t clustering 
requirement," JHEP 0212 (2002) 063 [arXiv:hep-ph/021 1426]. 

[16] A. Banfi and M. Dasgupta, "Problems in resumming interjet energy flows with k t clustering," Phys. 
Lett. B 628 (2005) 49 [arXiv:hep-ph/0508159]. 

[17] Y. Delenda, R. Appleby, M. Dasgupta and A. Banfi, "On QCD resummation with k t clustering," JHEP 
0612 (2006) 044 [arXiv:hep-ph/0610242]. 



6 



